BCS-BEC crossover in a two-dimensional Fermi gas 
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We investigate the crossover from Bardeen-Cooper-Schrieffer (BCS) superfluidity to Bose-Einstein 
condensation (BEC) in a two-dimensional Fermi gas at T — using the fixed-node diffusion Monte 
Carlo method. We calculate the equation of state and the gap parameter as a function of the 
interaction strength, observing large deviations compared to mean-field predictions. In the BEC 
regime our results show the important role of dimer-dimer and atom-dimer interaction effects that 
are completely neglected in the mean-field picture. Results on Tan's contact parameter associated 
with short-range physics are also reported along the BCS-BEC crossover. 
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The study of ultracold atomic Fermi gases has become 
an active and rich field of research Important ar- 
eas of investigation include the BCS-BEC crossover in 
a superfluid gas with resonantly enhanced interactions, 
the Chandrasckhar-Clogston instability of the superfluid 
state when spin polarization is increased, the possible 
onset of itinerant ferromagnctism in a gas with repulsive 
interactions Q and the realization of the Hubbard model 
for fcrmions loaded in optical lattices Q. 

Low dimensional configurations of degenerate Fermi 
gases have also been the object of experimental and the- 
oretical studies [H, In particular, a two-dimensional 
(2D) ultracold Fermi gas has been recently realized us- 
ing a highly anisotropic pancake-shaped potential, and 
the density profile of the cloud has been measured using 
in situ imaging [1]. On the theoretical side, the evolution 
from a superfluid state with large Cooper pairs to one 
with tight molecules in a 2D system of attractive fermions 
was first investigated by Miyake [H| and later by Randeria 
and coworkers [6| aiming to describe high-T c supercon- 
ductors. More recent studies address the problem of the 
superfluid transition 0, Q , of harmonic trapping 0] and 
of population and mass imbalance ■ These studies are 
in general based on perturbativc or mean-field (MF) ap- 
proaches that are suitable in the regime of weak coupling, 
but are bound to break down for stronger interactions. 

In this Letter we provide the first determination using 
quantum Monte Carlo methods of the equation of state at 
T = of a homogeneous 2D Fermi gas in the BCS-BEC 
crossover. We also obtain results for the pairing gap and 
the contact parameter as a function of the interaction 
strength. In the strong-coupling regime the emergence 
of interaction effects involving dimers produce large de- 
viations compared to MF predictions. A similar study 
carried out in 3D ll| has provided an important bench- 
mark against which experimental determination of the 
equation of state, using measurements of the dispersion 
of collective modes [H] or of in situ density profiles fl3| . 
have been successfully compared. Hopefully, our results 
will stimulate more experimental efforts towards the re- 



alization of a 2D Fermi gas in the strong-coupling regime 
by means, for example, of a Feshbach resonance to in- 
crease the interaction parameter jij. 

We consider a homogeneous two-component Fermi gas 
described by the Hamiltonian 



H 




(1) 



where m denotes the mass of the particles, ... and 
label, respectively, spin-up and spin-down par- 
ticles and Nf = N± = N/2, N being the total number 
of atoms. We model the interspecies interatomic inter- 
actions using an attractive square-well (SW) potential: 
V(r) = -V for r < R (V > 0), and V(r) = oth- 
erwise. In order to ensure that the mean interparticle 
distance is much larger than the range of the potential 
we use nR 2 = 10 -6 , where n is the gas number density, 
or equivalcntly kpR = 0.0025 in terms of the Fermi wave 
vector kp = \j2nn. We simulate a strictly 2D system 
and describe the low-energy collisions of the SW poten- 
tial in terms of the 2D scattering length <i2D defined as 
a2D = R e Jo( - K ^ Kjl( - K \ where J m(x) are Bessel func- 
tions of the first kind and k = ^Jv mR 2 /h 2 Q. The 
scattering length is non negative and diverges at k = 
and at the zeros of J\, corresponding to the appear- 
ance of new two-body bound states in the well. Close 
to these points the shallow dimers have size ei2D and 
their binding energy is given by Sb = — 4?i 2 /(ma 2 , D e 27 ), 
where 7 ~ 0.577 is Euler's constant The depen- 

dence of a2D on the depth Vq in the region where the 
well supports only one bound state is shown in the inset 
of Fig. [T] Two regions are clearly identified by com- 
paring a2D with the mean interparticle distance 1/Uf'- 
i) fcpa2D 3> 1 corresponds to the BCS regime where 
interactions are weak and dimers are large and weakly 
bound, ii) kpa,2B "C 1 corresponds to the BEC regime 
of tightly bound composite bosons. Compared to the 
3D case the BCS-BEC crossover in 2D exhibits impor- 
tant differences, a) For a purely attractive potential a 
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FIG. 1. (color online). Equation of state in the BCS-BEC 
crossover. Squares refer to the BCS and circles to the JS 
wave function. The solid (red) line is a fit to the data, the 
dotted (green) line is half of the molecular binding energy and 
the dashed (blue) line is the MF prediction. The horizontal 
dotted (black) line denotes the energy per particle Efq of 
the noninteracting gas. Inset: 2D scattering length a2D as a 
function of the depth Vo for a SW potential of radius R. The 
BCS and BEC regimes correspond, respectively, to fc_Fffl2D 3> 
1 and A;_Fa2D <C 1. 

two-body bound state exists for arbitrarily weak attrac- 
tions, b) The weak-coupling limit corresponds to a di- 
verging scattering length a2D- c) The 2D scattering am- 
plitude of particles colliding at low energy is given by 
f(k) = 27r/[log(2/fca 2D e 7 ) + i7r/2] 0. There is no range 
of values of a2D for which f(k) is independent of interac- 
tion (unitary limit), d) The mean- field coupling constant 
can be written as g = {2ttT% 2 /m)/ log(l/fci?a2D) with log- 
arithmic accuracy. Within the same accuracy, the region 
&f 0-2B ~ 1 identifies the strong-coupling crossover be- 
tween the BCS and the BEC regimes [see inset of Fig.Q]. 

Simulations are carried out in a square box of area 
L 2 — N/n with periodic boundary conditions, using the 
fixed-node diffusion Monte Carlo (FN-DMC) method. 
This numerical technique yields an upper bound for the 
ground-state energy of the gas, resulting from an ansatz 
for the nodal surface of the many-body wave function 
that is kept fixed during the calculation (see Ref. [l7[ for 
more details). The boundary condition is enforced u sing 



a trial function that we choose of the general form [18 1 
V't(R-) = $s(R)$a(R)- ®s is a positive function of the 
particle coordinates R = (ri, ...,rjy) and is symmetric 
in the exchange of particles with equal spin, while $^4 
satisfies the fermionic antisymmetry condition and de- 
termines the nodal surface of ipr- The symmetric part 
is chosen of the Jastrow form $s(R) = Yli v /ti( r "')> 
where two-body correlation functions of the intcrparticlc 
distance have been introduced for antiparallel spins. The 
<&A component is chosen as an antisymmetrizcd product 
$a(R) = A{(j>{r lv )(j)(r22')--4>{rN t N x )) of pairwise or- 



in the box, with integers ts summed up to the maxi- 
mum value of the fc-th shell accommodating N/2 par- 
ticles, and (3 is a variational parameter controlling the 
relative weight of the plane-wave sum to the spherical 
symmetric component <p s (r). Two important regimes are 
described by the above trial wave function: i) if /3 = 
and f s (f) — fb{f) is the two-body bound state of the po- 
tential V(r), V>t(R) describes a BCS state of dimers that 
is expected to be appropriate in the deep BEC regime; ii) 
if instead (p s = 0, the antisymmetric component in the 
trial function coincides with the product of the plane- 
wave Slater determinants for spi n-up and spin-down par- 
ticles, $a(R) = D^(N f )Di(Ni) (IsJ, and ip T is a typical 
Jastrow-Slater (JS) function of a normal Fermi liquid. 
This description is expected to hold in the BCS regime 
of a weakly interacting gas where the effect of pairing on 
the ground-state energy is negligible. The more general 
form of the trial wave function written above interpolates 
between these two regimes. 

In Figs. OH] and in Table U we report the FN-DMC 
results for the equation of state as a function of the in- 
teraction parameter in units of the energy per particle of 
the noninteracting gas E F q = h 2 kp/Am = Ef/2, where 
ep is the Fermi energy. Calculations are carried out us- 
ing i/jt of the BCS and JS form as described above. The 
BCS wave function (corresponding to /? = 0) provides 
a lower energy for values of the interaction parameter 
?7 = log(fci?a2D) 1, while the JS function (correspond- 
ing to j3 1) is more favorable for larger values of rj. 
The optimal parameter f3 in the BCS orbital has been 
found to be zero even in the region 77 ~ 1; finite val- 
ues of /3 have not given a significant improvement of 
the ground-state energy. The role of finite-size effects 
has been investigated by carrying out calculations with 
1 




bitals of the form 
k Q = 2Tr/L(£ ax x - 



r ) = PHk a <k F e a ' r + V>s{r). Here, 
t ay y) indicate the plane-wave states 
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FIG. 2. (color online). Equation of state in the BCS-BEC 
crossover with Eb/2 subtracted from E/N . Symbols are as in 
FiglU The solid (red) line is a fit to the data, the other dotted 
lines show the equation of state ([2]) of composite bosons and 
the perturbation expansion holding in the BCS regime (see 
text). The dashed (blue) line is the MF result. Inset: short- 
range behavior of g-fi for r/ = 2.15. 
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TABLE I. Energy per particle and molecular binding energy 
in the BEC-BCS crossover (energies are in units of -Efg). 



log(k F a 2 o) 


E/N 


£b/2 


E/N - e b /2 


-2.00 


-137.761(7) 


-137.832 


0.070(7) 


-1.50 


-50.593(4) 


-50.675 


0.082(4) 


-1.00 


-18.532(4) 


-18.637 


0.105(4) 


-0.50 


-6.714(4) 


-6.856 


0.142(4) 


0.00 


-2.318(2) 


-2.522 


0.204(2) 


0.25 


-1.283(12) 


-1.530 


0.247(12) 


0.50 


-0.638(10) 


-0.928 


0.290(10) 


0.75 


-0.201(12) 


-0.563 


0.361(12) 


1.44 


0.349(6) 


-0.143 


0.492(6) 


1.72 


0.459(16) 


-0.080 


0.539(16) 


2.15 


0.552(2) 


-0.034 


0.587(2) 


2.64 


0.634(4) 


-0.013 


0.647(4) 


3.34 


0.706(2) 


-0.003 


0.709(2) 


4.03 


0.755(4) 


0.000 


0.755(4) 


4.37 


0.775(1) 


0.000 


0.775(1) 


5.18 


0.821(7) 


0.000 


0.821(7) 



N = 26 and TV = 98. No significant change is seen when 
using the BCS trial function. In the case of the JS func- 
tion a large suppression of such effects is obtained by 
using the theory of Fermi liquids. The difference in the 
energy per particle between the finite-size and the in- 
finite system, in the interacting case, is assumed to be 
the same as in the noninteracting case, to lowest order 
in the effective mass (see |20j for details). The MF re- 
sult E/N = E FG + s b /2 (aB is shown in Figs. QUI for 
comparison. The inadequacy of the MF approach is best 
shown in Fig. [2 where the molecular contribution is sub- 
tracted from the energy per particle. This figure has to 
be compared to Fig. 5 of Ref. [l| , concerning the 3D case: 
effects beyond mean-field are much more pronounced in 
2D than in 3D. In the BEC regime the FN-DMC results 
are fitted with the equation of state of a gas of composite 
bosons corresponding to hard disks of diameter a d 
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2ttTl n d 
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log 7T - 



log(l/n d a 2 ) 
2^ + 1/2' 
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where m d = 2m is the mass of the dimer, while the num- 
ber of dimers, and correspondingly their density n dl is 
half of the total number of fermions A^ = N/2. The 
above expression includes beyond mean-field terms plj 
and allows for a precise determination of the dimer-dimer 
scattering length a d . We obtain a d = 0.55(4)a2D, in 
agreement with the four-body calculation in Ref. [7[. In 
the opposite BCS regime, where the contribution of the 
pairing gap can be neglected, the fermionic equation of 
state can be described in terms of an attractive nor- 
mal Fermi liquid (FL). Beyond logarithmic accuracy one 
has the second-order expansion in terms of r\ [22, HH 

jj = Efq (l — i + . From a best fit we find the re- 
sult A = 0.06(2) for the coefficient of the second-order 
term [H. 



In Fig. [3] we show the results for the pairing gap A gap 
in the strong-coupling regime. This quantity is defined 
from the difference of ground-state energy E(Nf, JVjJ 
of systems having one and two more (less) particles 
A gap = l/2[2E(N/2 ± 1, N/2) - E(N/2 ± 1, N/2 ± 1) - 
E(N/2,N/2)} 0. At the MF level 0, @ the pairing 
gap coincides with the result for the order parameter 
Ag a p = A = ^2£_f|£ 6 | if | Eb | < 2e_f, and is given by 
Agap = ef + |eb |/2 for larger values of \sb\- In the BEC 
regime the quantity A gap — |£f>|/2, shown in the inset 
of Fig. El displays the repulsive interaction effects be- 
tween unpaired fermionic atoms and bosonic dimers. In 
fact, the energy of the system with one extra spin- up 
particle can be written as the sum of the contribution 
@ of N/2 dimers and the Fermi-Bose interaction en- 
ergy E(N/2 + I, N/2) = E{N/2 1 N/2) + g BF n dl where 
9bf = 37r?i 2 /[mlog(l/nda^ d )] is the coupling constant 
fixed by the atom-dimer reduced mass 2m/ 3 and the ef- 
fective scattering length a ad - By using the definition of 
Agap and the value a d = 0.55a2D for the dimer-dimer 
scattering length in the energy functional ([2]), we find 
a ad = 1.7(l)o,2D from the fit shown in the inset of Fig. [3] 

Finally, we calculate the contact parameter C \2~h\ - 
[2^| defined from the short-range behavior of the an- 
tiparallel pair distribution function lim r _j.o 5tl( r ) = 
AC/kp log 2 (?'/a2D) [sec inset of Fig. [2]. The contact 
parameter is also related to the derivative of the equa- 
tion of state with respect to the interaction parameter 
C = {2Trm/h 2 )d{nE/N)/d{\ogk F a 2D ) The results 

are shown in Fig. |4j In the inset we show the quantity 
C — Co, where Co = (irm/h 2 )d(n£b) / 'c?(logfcFa2D) is the 
contribution to the contact C from the molecular state. 
The comparison between the two determinations of C is 
a stringent consistency check of the theoretical approach. 
We find a good agreement with Tan's relation, apart from 
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FIG. 3. (color online). Excitation gap in the BCS-BEC 
crossover. The solid (red) line is the MF result. Inset: ex- 
citation gap with |£f,|/2 subtracted from A gap . The dashed 
(blue) line is a fit using the energy functional of a Fermi-Bose 
mixture. 
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FIG. 4. Contact parameter in the BCS-BEC crossover. The 
solid line corresponds to the calculation from the derivative 
of the equation of state reported in FigQ] Inset: Contact 
parameter with the two-body contribution Co subtracted. 

the region rj ~ 1 where small deviations are visible, both 
with the JS and BCS-type wave function, showing the 
need of a better optimization of ipT- 

An important question relates to the relevance of these 
results for systems in harmonic traps. Two-dimensional 
configurations are realized if the transverse confinement 
is strong enough to reduce the kinematics to the xy- 
plane: Tiuj z Ef — TujJ±_-\/N , where we assumed isotropic 
trapping in the radial direction ui x = uj v = u±. In 
these conditions the effective 2D scattering length can 
be expressed in terms of a^o and the transverse har- 
monic oscillator length a z — y/h/muJz being given by 
a 2 D = Oz(2 A /7r/i?/e 7 ) exp(- > /7r/2q z /a 3D ), where B ~ 
0.905 0, [l6|. For small, negative values of the 3D scat- 
tering length the system is found in the BCS regime 
corresponding to an exponentially large 02D- The BEC 
regime is reached if the absolute value of <Z3D is increased 
such that |o3d| 3> a z j log(l/&Fa z ). An additional re- 
quirement concerns the dimer state, which is well de- 
scribed by the 2D expression only if \eb\ <C huj z [l6j], or 
equivalently a2D 3> &z ■ We believe that this latter condi- 
tion can be relaxed if, in the comparison with the results 
reported in this work, one considers quantities where the 
molecular contribution has been subtracted out. 
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work Programme. 
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